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This report is subdivided into three parts, each with its own page 


numbers, references, etc. ‘These parts are: 


18. Phase Shifter Develon:snt 
by Robert Seckelmann 


Tee Non-Linear Losses in Ferrites 


by Edward J. McKinney 


nea Propagation of TE Modes in Dielectrically Loaded Rectangular 
Waveguides 


ty Robert Seckelmanr 


Part I. contains data on the experimental work and general 


acknowledgements. 


Part II. and Part III. are theoretical studies. 


PART I, PHASE SHIFTER DEVELOPMENT 


ABSTRACT 


The phase shifters discussed are nonreciprocal. They contain circum- 
ferentially magnetized ferrite cylinders, which work at remanence. The 
differential phase shift is controlled by partial magnetization of the ferrite, 
achieved by controiled flux transfer. Three problems, common to all such 


phase shifters are dealt with: 


1. The ei.iciency of flux transfer and the possibility to add :lux bits. 


Experimental data and a theorv in agreement with these data are presented. 


2. The determination of r-f magnetic threshold fields above whith the 
losses of ferrites increase with field strength. The threshoid is cali:ulated 
from the measured threshold power of guides completely or in part filled 


with ferrite. 


3. The influence o* waveguide and ferrite dimensions on the frequency 
dependence of the differential phase shift, onthe proportionality betwen 
phase shift and transferred flux, and on the peak power handling capability 
is investigated experimentally. Slabline and waveguide type devices 


have been used. 
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I. PROJECT OUTLINE 


The object of the program reported here was to demonstrate the feasa- 
bility of: non-reciprocal multibit latched ferrite phase shifters for freq- 
uencies in the 4 - 12 GHz range, capable of handling peak powers of 10 KW 


and controlled by flux transfer techniques. 


Tentative goals were set in contract review conferences during the 


initial phase of the development as follows: 


Investigate slabline and waveguide type devices 

Switch in four bits (180°, 90°, 45°, 22.5°) 

Achieve at least 360°/dB 

Lrive with 500 wjoules or less per switching cycle 

Aliow for a switching cycle approximately 5 wsec. 

Understand the trade-offs between the various 

parameters and formulate design rules 

The phase shifter requirements could have been fulfilled by using 

four ferrite bits of different lengths along the axis of the device and 
switching them between their largest negative remanent magnetization and 
some positive remanent magnetization smaller than the largest one by means 
of controlled flux transfer to secure a phase shift independent of temp- 
apature ss) In review conferences, however, the investigators were 
encouraged to pursue their idea to achieve several bits of differential 
phase shift by using only one bit of ferrite alorg the axis of the device 
and to control the phase shift by controlling the state of magnetization 
of this ferrite. Consequentiy a device very close to an analogue phase 
shifter has been developed. This second approach requires a mvch tighter 
control of the field distribution in the device and imposes more res- 
trictions on the territe and guide shape than the first approach would have 


done. 


Flux transfer circuits have been analyzed theoretically and the theory 
has been found to be in good agreement with experimental data. Circuits 


have been designed to add phaseshift bits by adding flux bits. 
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The phase: shifters and the syst. us using them could be considerably 
simpler, if it were possible to eliminate some of tlhe switching wires. 
This should be possible by supplying to the phase shifter voltage pulses 
of constant height (voltage) and variable length (time). Ideas ror such 
circuits - temporarily called “analogue drivers" - have been conceived 
during the last phase of the contract, but have not yet been tested. An 
electronic driver of very small volume and steered by a computer could be 
incorporated ir the phase shifter. Whereas now such a computer has to 
select - in the case of a four bit phase shitter - any of 16 bit combinations, 
it, would then select any of 16 times, during which it would activate the 
driver. All phase shifters studied use round ferrite cylinders. The 
devices have too complex a geometry to be directly dealt with mathematically. 
Rectangular waveguides with dielectric slabs parallel to the narrow guide 
walls serve as more or less crude analytical models. ‘The propagation of 
TE modes in these loaded guides has been studied in detail. Nume>ical 
results for a large range of parameters are given in Part III of the 
report. They will also be published as a G.E. TIS Report. The computer 
program <o obtain these data was financed in part by this contract. The 
Slabline type phase shifter performance could be explained nicely in terms 
of TE modes in dielectrically loaded rectangular waveguides. The computer 
results were also used to evaluate tests to determine the threshold field 
strength for the onset of non-linear losses in a rectangular waveguide 


containing a ferrite slab. 


The differential phase shift versus flux dependence led to a model for 
magnetization of ferrite toroids, which in turn led to @ reciprocal TEM 
phase shifter controlled by partial magnetization. This reciprocal phase 


shifter has been tested in one of the high power test structures. 


The development of high-power ferrites is covered by other contracts. 


Here the influence of the device geometry on the on-set of non-linear losses 


was considered. ‘he requirements were met, and results encouraging the develop- 


ment of phase shifters capable to handle higher peak powers, qey 50-100 KW 
were obtained. An effort has been made to understand and predict the onset 


of non-linear losses in ferrites of arbitrary shape. A theoretical study is 
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given as part II of this report. This urcatise includes the results 
obtained by other workers with more restricted models, At the moment, 
however, it cannot be used to its full extent because not all parameters 


could be determined in experiment. 


Figures of merit of more than kOO°/ dp have been obtained throughout 
the development. Nearly octave bandwidths with only a few percent variation 
in differential phase shift have been obtained in waveguide type phase 
shifters. The nighest peak power handling capability, the broadest band- 
width and proportionally between differentia: phase shift and flux 
transferred have, however, not been obtained simultaneously. Trends have 
been found for the optimization of each characteristic, but the trade-offs 
are not fully understood. Nearly all parameters are nutually dependent 
(e.g. ferrite wall thickness, diameter. wavezuide height and width, 


filling factor). 


The project was started with three parallel efforts: 
to investigate flux transfer circuits 
to study the onse. of non-linear losses in ferrites 


to develop phase shifters at low powers. 


™T the second half of the project the flux transfer investigation was 
replaced by an effort 


to develop criver circuits 


Finally, the phase shifters were tested with high peak power pulses. The 


different efforts will now be described in the order given above. 


II. CONTROL OF FLUX-TRANSFER 


A flux transfer circuit is shown in Figure 1. R is a variable resistor in the 


circuit. Cp © Cy, ..., Cp are n “driver” cores (index D) of various sizes, which are 


switched from negative to positive saturation magnetization or the corresponding 
remanent magnetization. CC is the "microwave" core (no index), which produces the 


differential phase shift. For both, the driver and the microwave core the simple BH- 
loops of Figure 1 are used to analyze the transfer circuit. B is the flux density and H 


is the magnetic field strength. 


In an unmagnetized ferrite core the magnetic dipoles are randomly oriented. 


exist domains for any orientation. Ina core saturated by sending a etrong current 


pulse through the core the ferrite is circumferentiaily magnetized. 


Sending now an 
increas!ag current in opposite direction through the core, one reverses gradually the 

sense of circumferential magnetization. Two models to describe the magnetizing process 
may be used, Ina "many shells" model the ferrit~ may be looked upon as consisting of 


many successive, very thin, laminar layers. A magnetizing current will either not 
change the state or magnetization of a particular shell or drive it to saturation 
magnetization. According to this model, an increasing d.c, current would at first 


"awitch” the inner ferrite rings, then the outer ones. An a.c, current with decreasing 
anplitudes would produce a sequence of circumferentially magnetized ferrite shells with 

alternating directions of magnetization. In a "one shell" model the magnetization would 
be homogeneous throughout the ferrite wall and can vary from complete spin alignment in 


me direction to random spin orientation to complete alignment in the other circun- 


ferentia’ direction, depending on the driving conditions. Ana.c. current pulse with 
decreasing amplitudes sent through a "switching" wire through the axia of the ferrite 
It is felt that 
often a super position of: a "many shells" and a "one sheii” »nodel will offer the best 


core wouid demagnetize the ferrite, i.e. randomly orient the spins. 


explanation of effects. Assuming for the ferrite a simplified’ BH loop as shown in 


Figure 1, the shell thickness would te given approximately by ry'/ro' = 
ry',.To' are the radii of the layer switched (with ry) Sr)' s Po! S ro and rj, rp the 


diner and outer core radii), 


Dynamic permeabilities pu = 4B/dH are defined by — 
pO) = B/Uy — HDa Hy * Byy/Clgg = Hye) 


(H) . (H) 
(B, ~ B,)/Hy» Hp (Bh, 3 By.) /Hyg 
and dynamic inductances with the core height h by 
L = df/aI = yh ln (ry'/P,')/en 


H, where 


(1a) 


(1b) 


(2) 


A' = alr, ! - r,') is the area of a particular core, in which the flux it Changed by 


= «Al i] zs 
4p = AB A, with Ai = A. 
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FLUX TRANSFER SWITCHING TECHNIQUE 
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Before a switching cycle begins, a reset pulse drives both the driver and the 
microwave core to saturation and subsequently to the corresponding remanent state. Let 
this be the negative state for both cores, - B, and -Bp., respectively. The driver core 


HET 


1s now driven to Bog and then falls back to Bp, = Bos - Bp: the remanent state, It 
feeds the transfer circuit with an approximately rectangular voltage pulse of height 
(voltage * V, and length (time) T, followed by a pulse with a relatively small negative 
peak voltage. T is proportional to the length of the driver core. The avaiiable flux is 


Of = 2B, A= Af) - bf), Og, = | Vi(t)at po Vio (3) : 
fo) 
Let 8 << 4g, (with holding current applied or uy es ()), so that 
' 
‘fp om Of, Cs) ; 
The flux transferréd is fo) 1 4 
4g. = J Va (tjdt = | Vi(t)dt - 6¢ (5) 2 
t 6 2 2 = 
6f is a relatively small contribution caused by the falling back of the microwave core a. 
from a value B, to By - 6By. With pp nw O,the time constant for this fall back : 
process is given by L(H) pp, Let the microwave core be driven to a value B< Be: The _ 
transferred flux is for a lossy circuit: 2 
a F 
A. sy | vi(t)dt - 4(t) Rat (6) a 
, i 1 = 
3 
Eo q 
ap, + of = | Vi (t)dt - J s(t )Rae 2 | Visite) ae (7) i 
so that : 
i(t) = (1 - exp(- tR/L(t))) V,/R ow tV,/L(t) with t s T << L(t)/R = 
Initially, during 0 st s Ty << T, the inductivity of the micrcwave core is LtH) The = 
current in the transfer circuit rises fast as 
I(t) = (1 - exp(-tr/i‘4))) V,/R ~ ev, /ult) witht <3 <<7 (Ba ) ‘ 
x 
At t = T, the inductivity becomes LV) >> pH) &ud the transfer current increases only = 
slowly a = 
Ut) = 17) (1 = exp(-(t = 2,) w/t”) (W, = 1(2,)R)/R $ 
Tt + (Ce - 2); I(T, )R)/L I(T, ) (8p) = 


with tT, sts p<) yp 


= 
= 
3 
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A HP 214A pulse generator was used to saturate the driver core. 


Figure 2 shows the output voltage of the unloaded source, the output 
voltage of the driver-core-loaded wource, the output current of the driver- 
ecre-loaded source, the voltage across the secondary side of the driver ° 
core (open transfer circuit), the voltage across the secondary side of the 
driver core (with microwave core applied). One recognizes the fast increase 
in current (from the source) due to the small inductivity between - Bo. and 
Boye then the slow increase of current due to the large inductivity between 
* Boe and Bs an. the fast increase between Bos and B > Boe Figure 3 shows 
the current in the transfer circuit (R = 0) and the voltage across the micro- 
wave core. The decaying current trail at the end of the current pulse 
produces a negative voltage across the microwave core, as indicated in 
Figure 1. in Figure 2 it cannot be seen. It follcws the positive voltage 
pulse and has a very low peak value. Figure 4 shows the voltage vy across 
the driver core for several] driver core sizes (with holding current). The 
driver core supplies an approximately rectangular voltage pulse of height 
(voltage) vs, and a width (time interval) T. In Figure 4 vy is approximately 
comstant and T proportionsl to the size of the driver core. It is not 
necessary to produce the time voltage area vit via a driver ferrite. It 
can also be done electronically. As indicated in Figure 1, the flux loss 


via ohmic losses during the current rise time T, is very small compared 


ag 


witn the loss during T, st S$ T. With this approximation and T - T, eT, 


Ny 
the transfer efficiency becomes 


T= (A6,/08,) a 1 - (69/a9,) - (rr, /i'”?) (9a) 


(Vv) 


As long as tne microwave core is not saturated L* ‘wz bp»: In the "many 
shells" model this means that only a certain fraction of the tota] number 
of laminar layers are switched; in the "cne shell" model ee & function 


of an effective permeability, which increeses with driving current. 


In Figure 5 tue driver and microwave core sizes have bee varied by 


using various numbers (np and n) of equal ferrite rings for both, the 
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values ‘that the transt rv effleiency decceases somewhat. with Tnocrcasint step 
Pmbers Tae ey ie Pel Une Tike per ale is ploled: vyerauseti Ee total 
flux transferred. The results of Figure 6 indicate ‘that the 
vertical Drier “tbe Bit Poop Ss nete pare ied 2a the fauxie,, As secn 


1} 


later, Wis @Preeo can be conpensated Ean LHerease Wir aseushiit tor 

&@ constant fiux transferred when stepping up the Hli-ioop. Ina first 

approximation hewever, the transfer exfi-ciency can be jooked upon as 
onstant until the micrcwave core is nearly saturated. The transferred 
flux is then by, - 64. If two bits Ay and 4. are switched one after the 


other, the total flux transferred is 


(2. + 2) (2 
Bee bee * af, = af, + Of, - 266 (11) 


If the two bits are switched at a time, the transferred flux is 
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TRANSFER EFFICIENCY WHEN "STEPPING UP" A BH-LOOP 
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The difference can be overcome by adding to the nr different bits a 
ccapensation flux bit by. = 64. Whenever bits are switched, iIndividuatly 
or several at a time, the extra bit is switched at the same time. The 


transferred flux is then 


09°) = (af, + 06.) - 66 = ag (13) 


Op, %7)- (ap, + a6.) - Of + (ad, + ag.) = Of = ag, + af, (24) 


cl ail aa~! ea guia hehe, 


af?) . (ag, + af, + af.) - Of = Of, + f, (15) 


The transfer circuits of Figure 1 and Figure 8 contain the microwave core, 
the active driver core to achieve a certain flux transfer,and some idle 
driver cores. The influence of the idle core on the transfer efficiency 
has been found to be negligible. ina setup as in Figure 1 different 
numbers of equal ferrite toroids have been used for microwave, driver and 
idle cores to establish this result. The idie cores are, of course, reset 
to negative saturation just as the active driver before flux is transferred. 
Thus, their inductivity is close to zero. The results are sumnarized 

in Table 1. The flux transferred to both the microwave and the idle 

cores ..as been measured. The low efficiency for o= 20% is due to the 
fact that rather small cores have been used, so that ohmic losses in the 


transfer circuit represented a considerable fraction of the total losses. 


Active [eae Idle Flux Microw. Flux 
cores avallable cores Transf. cores transf. oO N 
5.32 uVs 0 0 uVs 5 4.46 pVs 40% 83.7% 
ye Suen 3 “On02. 5, . HS HOU 36% 
J: Pig e\cta a @) @) 5 1.66 " 20% 62.5% 
} 2G! 4 <G) je 5 1.66 "20% 62.5% 
TABLE 1. Influence of Idle Cores on Transfer E1 ficiency 
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TRANSFER CIRCUIT FOR ADDING FLUX BITS 


FIGURE 8 
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j III. NON-LINEAR LOSSES IN F’RRITES 
{ 


Non-linear losses occur when the rf - magnetic field in the wave- 
guide surpasses a threshcld value. Power is then dissipated into the spin 


wave spectrum. The dominant mechanism is subharmonic generation, where 


spin-waves of half’ the frequency of the drivin. rf-field are excited. This 


j is a "first-order" effect because the spin wave anplitude is prcportional 
to the rf-field amplit:de. Via a second mechanism spin-waves can be ex- 


cited at the frequency of the driving rf field. This effect is of "second 


order" because the amplitude of these spin waves is proportional to the 


square of the rf-field amplitude. Whenever the subharmoni: generation takes 


place it occurs at such low rf-field amplitudes that the second order 


asia 3 


nes derived threshold field strengths for ferrites of arbitrary shape at 


remanence. This treatise is given fully as Part II of this report. It 


lala 


includes the :esults of other workers. McKinney has to normalize the 
threshold field strength to a spinwave line width or relaxation time. 
Were this line width known, it would be possible to predict the threshold 
field sirengtn tor ferrites as a function of frequency, and for devices 


where the field distribution is known one could also predict the threshold 


at one frequency only, even less to determine whether the relaxation time 
is a function of frequency, ferrite shape, driving field strength or other 
parameters. 


By testing ferrites of equal size and shape in equal test devices at 


series were per*ormed: In the first one a coaxial line was completely 


filled with ferrite. The diameters were chosen so that only TEM modes 
transformers at both ends. The VSWR was below 1.5 in all cases, correspendi 


to a reflection loss < 0.17 dB or a power reflection of iess than 4%. In 


the second series a ferrite slab was placed in the center plane of 1 rectan 


ing 


erst pig grag ily ils iis ii 
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effect may be neglected. Assuming only subharmonic gmeration, E. McKinney 


power. So far it was not possible to measure the relaxation time accurately 


a5 GHz, relative power handling capabilities could ve determined. Two test 


should propagate. This test device was matched into 50 ohm lines via step » 


ne 


gular 


> mane 


waveguide parallel to the narrow well and extending over the full height 

of the guide. The filling factor was five percent. The ferritc was tapered 
at both ends to secure a VSWR < 1.5. The rf magnetic field amplitude was 
calculated under the assumption that the relative ferrite permeability Hp, = 1 
and the relative ferrite dielectric constant Ep = HU, 25.16, 46 given vy the 
manufacturer. In the coaxial line the strongest rf-magnetic field occurs 


at vhe inner conductor. With 


Tecate rt 


P = 21° = 2m yp? # (ae) 


where P = peak power, I = current in inner conductor, D = diameter of inner 


conductor, Z = line impedance, H = rf-magnetic field. One obtains with 
2 


nm f& 10 
ee (17) 
A/cm cm Z/ohm fi 
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In the rectangular waveguide (WR 137) av the test frequency of 6.5 GHz the 
field distribution cucside the ferrite slab follows a nyperbolic function 


and one obtains with the dimensions sketched below 


W 


= (E 
P= (E. x H) dA = 2b J (Ex H, Jax 


(18) 
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wu me 
2wb se i ee ap + J eine (QP + 6) a? 
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(x) 
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where 


k a longitudinal propagation constant 
P/w 2 transverse propagation constant in empty part of guide 
Q/w © transverse propagation in dielectric 
1/D 2 relative amplitude in empty part of guide 

1 2 relative amplitude in dielectric 


6 = aw 
oe eg An 


The full theory is given as Part III of this report. For the given frequency 
and waveguide (Qw = 1.372 inches, b = .622 inches) one obtains finally 


any a (19) 
A/em tt V kW 


where N is a function of wne ferrite dielectric constant 


bit 1 8 ; ie 14 16 


1925. 15055 1.108 1.136 1.162 1.180 


The test results are summarized in Figure 9. The coaxial line test leads 
to considerably higher values for the threshold field than the wave guide 
test does. The reason is not yet understood. The spin wave spectra excited 
in both test series may not be equal. ther possible explanations are for 


the coaxial line test: 


The impedance is other than calculated with Ee 


Several modes exist (partly due to air gaps 
between ferrite and conductors) 


The rf-magnetic field strength varies con- 
siderably over the ferrite wall thickness. 
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Ferrites at High Peak Powers 
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The average threshold power in the waveguide test could be determined 
accurately, because the non~linear lose/power characteristic is a straight 
line in a double logarithmic scale as shown in Figure 10. If the ferrite 
wou d not warm up, this line would end abruptly at a maximim power-independent 
loss. Such a maximum loss nes been maintained over a large power range 
when the ferrite was kept cool. Figure 10 shows a peak in lcss due to the 
warming up of the ferrite and the subsequent change in saturation magnetiza- 
tion. The onset of non-linear losses in the ccaxial line tests was not as 
well defined. The losses inc.‘eased to a small. value over e large power 
range before they followed the straight line. Besides the results were not 


as reproducible as in the waveguide tests. 


Figure 11 shows the input pulse an@ the deformed output puise at 
various power levels. The reduction in loss at high powers shows up clearly. 
To determine the threshold power it has beer assumed that the input pulse 
is rectangular and 0.42 psec wide. Assuming an exponential growth cof whe 
spin waves it should be possibie to find the relaxation time from the decay 


of the high peak of the output pulse. 


With the poor input pulse shape, however, this does not yield accurate 
results. The spin wave line width appears to be between 0.5 to 2 A/cm. 


ile SMA 


The resonance frequency has been taken from the manufacturers data or 
from the onset of losses at lower frequencies in a coaxial line filled with 


ferrite, as illustrated in Figure le, or from measuring the flux available 


sj 
from a toroid. In any case the accuracy may be + 5%. Plotting the threshold 3 
fields versus w /», as done in Figure 9 and assuming a linewidth of 1 A/em, a 
one finds that the waveguide test data come closer to the theoretical results 
than the coaxial line test data do, 
In any case the power handling capability increases as the saturation ES 


ma,netization decreases. 
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IV. RECIPROCAL TEM MODE PHASE SHIFTER 


The "one shell" and the “many shells" model for vhe magnetization of 


ita 


ferrite voroids - mentioned under "Control of Flux Transfer" - leads to 

an interesting consequence. Tne "many shells" model forbids, while the 
“one shell" model allows @ reciprocal, latched ferrite phase shifter using 
the TEM modes in a coaxial line. According to the "many shells" model 

an alternating current with decreasing amplitudes applied to the inner con- 
ductor of the coaxial line would produce a sequence of circumferentially 
magnetized ferrite shells with alternating directions of magnetization. 

In any such sequence no microwave interaction between the ferrite and rf 
magnetic field of a TEM mode would occur. The effective relative permea- 
bility for this mode wculd always be one. Ina "one shell" model the mag- 
netization would be homogeneous throughout the ferrite wall thickness and 
could vary from complete alignment in one direction to random spin orien- 
tation tc complete alignment in the other circumferential direction, 
depending on the driving conditions. This model allows a reciprocal latched 
coaxial phase shifter. The -lectrical length of the device depends on the 
degree of dipole alignment in the ferrite, not on the circumferential 
direction. In a coaxial line filled with ferrite (TT1-109, L = 4 inches, 
D, = .221 inches, D, = .063 inches), a differential phase shift has been 
measured between a) the demagnetized state, achieved by sending an a.c. 
pulse with decreasing amplitude through the center conductor and, b) 

the state of partially circumferential magnetization (either direction) 


achieved by sending a d.c. pulse through the inner conductor. 


The reciprocal phase shift could be controlled by the amplitude of the 
d.c. current pulse or by controlled flux transfer. At about 5 GHz the 
reciprocal phase shift could be stepped up to approximately 20, 40, 60, 
and 80 degrees. The coaxial line used in this experiment was built for 
ferrite tests at high peak powers at 6.5 GHz. The phase shift experiment 
was only a feasibiiity demonstration. The demagnetization was done at 
60 Hz. The switching properties of nonreciprocal phase shifters - investi- 
gated under the present contract— indicate that demagnetization at a few 


kHz should te possible. 


raat 


Stepping up the B-H loop from one state of circumferential magneti- 
zation to the demagnetized state to the opposite state of magnetization gave 
a@ smaller maximum reciprocal phase shift than starting from the a.c. 


demangetized state. 
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V. DIFFERENTIAL PHASE SHIFTERS 


Three types of nonreciprocal phase shifters are shown in Figure 13. 
All of them use as the nonreciprocal element ferrite toroids, which are 
circumferentially magnetized and work at remanence. In all of them the 
phase shift is controlled by partial magnetization, achieved through 


control’.ed flux transfer. 


The "nurmal" helix type phase shifter has been developed and described 
Sy HAs Hair, (2) The "inverted" helix type phase shifter has been found 


less effective (2) than the normal one. 


The slabline type phase shifter and the waveguide type phase shifter 
with TJerrite toroids were uszd in the experiments described further down. 


Rectangular ferrite tubes were used in the first waveguide type phase 
f9\ 
shifter by levey and Silber’! 
(5) 


and the improved versions by Blevins, 


and Taft and Sweeney. 


The direction of magnetization in the phase shifters considered is in 
a plane transverse to the direction of propagation. Differential phase 
shift can be obtained when circularly polarized rf fields exist in a plane 
transverse to the direction of magnetization. Hybrid modes, which do not 
occur in empty guides, and TE modes, which have counterparts in empty 
guides, fulfill this condition. The tensor permeability of the ferrite 
makes the exact computation of the propagation constants in “hese guides so 


difficult, that so far no exact solution exists. An analysis has to be 


based on more or less crude models. In simple configurations exact solutions 


may be found for aielectrically loaded guides. Perturbation calculations 


(6) 


configuratiorB investigatec, the dielectrically loaded waveguide model 


ean then take care of thin magnetized ferrites. But even in the complex 


anaiyzed in Part III of this report offers food explanations. 


In the phase shifters considered here it has been attempted to achieve 
a proportionality between tue dif*erenti. | phase shift and the flux trans- 
ferred to the phase shifter. The ferrise was driven between its maximum 
negative remanent state of magnetization and various other remanent states 
between this one and the maximum positive remanence. All tests were per- 


formed with Trens-Tech ferrites. 
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Val Waveguide Type Phase Shifters 


All the phase shifters investigated use round ferrite cylinders. 
Previous work dene in this laboratory indicated, that they should perform 
at least as well - as far as figure of merit and bandwidth are ccncerned - 


as phase shifters using rectangular ferrite cylinders. 


The first experiments were done witn a TT1-105 ferrite toroid 
(M, es a7 us /om”) with .662" outer diameter in a standard WR 137 wave- 
guide. With tapered slabs the guide width was reduced. The final step was 
to go to a circular guide of the ferrite diameter. Tne device is shown in 
Figure 14. Other parameters varied were the ferrite wall thickness and the 
core dielectric constant. The results ure shown in Figures 15, 16, es 
(phase shift versus frequency) and Figure 18 (phase shift wnen stepping up 
the B-H loop). The results show that thin wall ferrites with proper di- 
electric cores can produce phase shift nearly as large as achieved with thick 
wall ferrites and that the phase shift versus frequency curve is the flatter 
the wider the guide. The "phase shift per step" characteristic for stepping 
up & BH-loop could be improved when the guide cross-section was changed as 
shown in Figure 19. The guide height is reduced +o less than the ferrite 
toroid diameter and the toroid rests in groves in the top and bottom wall 


of the guide. The threshold power with this device was /6& KW. 


Since the threshold field strengths fcr ferrites generally increases 
with decreasing magnetization the ferrite TT1-1400 (M, = 10 pVs/cn”) was 
chosen for further experiments. Figure 20 shows first results with this 
ferrite. Reducing the ferrite diameter reduced the slope of the phase 
shift versus frequency characteristic. Essentially, constant differential 
phase shift over a bandwidth of 2.5 GHz, centered at 6.2 GHz was achieved 
with a filling factor of about 31%. With increasing filling factor the slope 
becomes negative, as seen in the illustrations of this report; with de- 
creasing filling factcr the slope becomes seen as shown for guides with 


(5 


rectanguiar ferrite toroids by Taft and Sweeney. The optimun filling 
factor for constant phase shift over a broad band seems to be in the neighbor- 


hood of 35%. In Figure 22 it is seen that the guide height has no great 
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PHASESHIFT IN WAVEGUIDE WITH REDUCED WIDTH (Ws 782") 
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influence on the phase shilt. Imncreasimg the wall thickness by 80% 
increased the phase shift by 40%. The threshold power was independent —{ 


guide height and ferrite wall thickness. 


The flattest phase shift versus transferred flux characteristic at 
6.5 GHz was achieved with a filling factor of about 36%, as show in the 
center row of Figure 21. The pertormance otf this phase shifter is illus- 
trated in more detail in Figures 23, 2, and 25. ‘The phase shift/flux 
curve (Fig. 24) indicates that a "one driver bit" phase shifter could be 
built, if enough time is available to step up the RH-loop. The driver bit 
would be charged and discharged to the phase shitter as often as necessary 
to achieve a certain phase shift. The high power rroperties shown in 
Figure 25 are tyrical for all the phase shifters investigated. The ferrite 
was driven tetween its maximum remanent states. ‘The differential phase 
shift stayed constant up to a power 1 to 3 dB above the threshold power. 
Then it decreased. This decrease is mainly due to the ~hange in insertion 
iengtn in one state of magnetization, the length in the other state stays 


approximately constant. 


If the ferrite warms up even without non-linear losses, the insertion 
length and differential phase may vary even then. Methods to prevent this 
change are a) the uge of temperatu-e compensated garnets, b) an effec- 
tive cooling method, and c) to drive the ferrite between its maximum 
negative remanent magnetization and a positive magnetization smaller than 
the maximum. This last method has been demonstrated first by Hair, (1) 

It is inherent in the flux transfer controlled phase shifters investigated 


here, 


While this report was being written a final test was made with four 
inches of TT1-900 ferrite (M, ~~ 7 uVe/om"). The ferrite was selected 
because it was found to possess a very high threshold field strength 
(see Section III). The maximum differential phase shift was 104°, The 
bandwidth and flux contro]. of phase shift were comparable to the result 
achieved with TT1-1400 with the same device and ferrite dimensions. The 
onset of non-linear losses occurred at above 100 KW. At 150 KW the non- 
linear losses were < 0.5 dB. 
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High Power Performance of Differential Phase Shifter at 6.5 GHz 


FIGURE 25 
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Vise Slab Line Type Phase Shifters 


A typical differential phase shift versus frequency characteristic of 
a slabline type phase shifter is shown in Figure 26. The phase shift 
increas*s gradually with frequency below / GHz, shows a high peak at 
7.5 GHz and ther: declines suddenly. The phase shift per step of this 
phase shifter (with a reduced ferrite length) at 6.4 GHz shows a peak at 
the eighth step. 


With some imagination one may use a rectangular waveguide with di- 
electric slabs parallel to the narrow walls as a model to analyze the 
fields in the slab line. This is done in Figure 27. The TE) 5 mode for 
the centrally loaced guide is the first to appear. With increasing freq- 
uency the plane of circular polarization moves from the center of the 
dielectric towards its outside. In the slabline this is the region of 
magnetization. A "positive" phase shift increasing with frequency is to be 
expected. A maximum differential phase shift per step will occur when the 
plane with circular polarization is magnetized. At high frequencies this is 


4 


the boundary pla:re" tween dielectric and the em‘y part of the guide, 
coriresponding to the outer ferrite diameter. At these frequencies, however, 
a TED) mode corresponding to the side wall loaded rectangular guide is 
possible and at little higher frequencies the TEQ5 mode for either the 
centrally or sidewall lcaded guide may also exist. Both these modes 

Froduce close t> their respective cutoff-freqiencies "negative" phase shift. 
At nigh... frequencies negatively as well as positively circulariy pclarized 
fields exist in the dielectric. The slabline type phase shifter responds in 
this frequ ncy range with a rathe: erratic phase shift versus frequency 
characteristic. The phase shifter should, therefore, work in the freq- 
uency range below tne cutoff-frequency corresponding to the sidewi.1* 


loaded guide. 


Figures 28 ard 29 show that all ti.e slab line type phase shifters 
with parallel side wails have the frequency characteristic iescribed above. 
Increasing, in Figure 29, the wa]l thickness by 80% incr :ased the phase 


shitt by 504 at 6 GHz. ‘The threshold power ccuubled from 5 to 10 KW. 
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The lower value achieved with the thinner wail is somewhat doubtful because 
moding problems occur at the test frequency. In some other slabline type 


phase shifters with ferrites with D. = .5", D, = .3" the threshold powers 


6 
achieved at 6.5 GHz were 
6 KW peak with TT1-105 in a shielded slabline 
9 KW peak with TT1-1C5 in an open slabline 


15 KW peak with TT1-1400 in an open slabline 


cr eh 


Changing the ferrite and device dimensions led to the phase shifter of 
Figure 30. The figure of merit is above 400°/aB. Half the ferrite is in 
contact with the outer conductor. This facili+ates heat transfer problems 


at high average and peak powers. Figure 31 ows the phase shift per 


step for this device for two different bit sizes and their combinations 
(switched one after the other and both together) with and without a 
correction bit. Though the device performance is not yet ideal, the 


influence of the correction bit is striking. At 5.7 GHz four equal phase 
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shift bits and a correction bit have been added with a transfer circuit 


as in Figure 8. The results are: 


Bit dp/deg(measured) d&p/deg(nominal) 
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Adding Phase Shift Bits by Controlled Flux Transfer (f = 5.6 GHz) 


FIGURE 31 
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Figure 32 shows the results achieved with four different bit sizes and the 
correction bit. The bits are somewhat too small for a 360° phase shifter. 
The best approximation is given by the 350° line. With respect to this 

line the error was within + 6° and within + 3%, whatever is smaller. Better 
results have been achieved but wer+ not reproducible. For good results 

the switching wire has to be perfectly parallel to the center conductor 

and center conductor and switching wire have to be perfectly centered. This 
phase: shifter could not be tested at high power, because the available 


source works at 6.5 GHz, i.e. beyond the phase shifter range. 
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To enlarge the bandwidth of the slab line type phase shifter a thin 


slab was inserted into the slab line as shown in Figure 33 to suppress the 


alia | 


mode corresponding to the TE,) 


guide. The frequency range has well been extended to higher frequencies. 


mode in the sicewall loaded rectan,ular 


The small remaining erratic phase shift response is ascribed to imperfect 
adjustment of the switching wire and to the mode corresponding to the TEQo 
mode in the rectanghlar eielectrically loaded guide. Non-linear losses 

occurred at about 2.5 KW and voltage breakdown at about 35 KW. The figure 


of merit is > 300 deg/dB. 
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32. Slabline Type Phase Shifter with Flux Controlled Phase Bits 
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VE General Remarks 


The differential iine Jengths is 


AL = (tp/rad) . (A/2n) 


AL/om = (1/12) . (dp/deg)/(£/GHz) 


Phase shifters with positive slope [a(dp)/df] may thus be useful rs delay 
lines with frequency independent differentiai line length. All slabline 
type phase shifters investigated have this positive slope. While the 
differential phase shift of the phase shifter in Figure 30 increases by 
48% between 4.7 and 5.€ GHz, the differential line length increases only by 
15%. The waveguide type phase shifter exhibits no slope with a filling 
factor of about 1/3, positive slcpe requires a reduced filling factor. A 
very broad band with constant phase shift has been achieved with this phase 
shifter. For low pover applications the small size of the slab line type 
phase shifter is an advantage, for high power application voltage break- 
down ma: occur before the ferrite goes unstable because of the small 


dimensions. 


How well a ferrite is used in a phase shirter at a certain frequency 
may be described in terms of phase shift per flux. Typical for the devices 
investigated is 10° /uVs. This value may be constant over a broad band. 

In helical phase shifters one may obtain 30°/uVs in the narr w band with 


maximum phase shitt. 


7) 
Other workers, ufing rectangular ferrite cylinders in waveguide 


type phase shifters, found that rounding the edges of the ferrite cylinders 
increased the differential phase shift and the figure of merit. Under 


the present contract the "rounding effect" was carried to the extreme. 
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VI. DRIVING CIRCUITS 


Figure 34 shows the basic circuit used in the drivers. It produces 
an eleven ampere current pulse 2.5 sec. .ong with 0.5 psec. rise time and 


0.8 usec. fall time as sketched below. The power supply voltage was made 


Trt \ 
vs a 


re eR = t/s 


as large as possibile under the constraints of the 2N1908 transistor. 


Initially the reset driver was arranged as in Figure i. ‘This arrange- 
ment did not work well because one of the two parallel inductivities 
(microwave and driver core) wag: saturated first and then shorted the other 
one, so that the other core was not completely reset. Therefore, an 
arrangement , as indicated in Figure 8, was chosen for the reset circuit. 

A block diagram for the complete driver circuit is given in Figure 35. 


The phase shifter of Figure 30 has a ferrite volume VAs 14 cm”, a 


coercive field strangth H, 1A/cm, and a saturation magnetization 

M, # 10 uVs/em@. The required energy per switching cycle is 

E equ oad 4H MoV 99560 pjoules. The drivers have an efficiency of at 
least 65%. Thus the total switching energy per cycle is 840 pjoules or 
less. It is believed that :for all the phase shifters investigated one 


needs 0.8 tc 1.5 pjoules/degree. 
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EXPERIMENTAL FOUR BIT DRIVER CIRCUIT 


FIGURE 35 
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PART II, NONLINEAR LOSSES IN FERRITES 


ABSTRAC 


First and second crder solutions for the critical field at which 
the magnetization in a domain becomes unstable are derived. The domain 
is assumed to be immersed in @ larger sample of ferrite which has a spacial 
average magnetization equal tc the remnant value. it is assumed that the 
magnetizaticn in the domain of interest is uniform and in the direction 
of the applied field since this case yields the smallest critical field. 
The first order solution corresponds to the zero order spin mode for the 
domain (all spins precessing in unison). The second order solution 
includes spin waves which exist as traveling waves within the domain. A 
third order sclution would include interactions between spin wave modes, 
particularly the zero order mode and higher order modes. The spin wave 
which is excited to unstability by the driving field is shown to have 
half the frequency of the driving field. It is shown that the critical 


field increases as the saturation magnetization of the material decreases. 
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NON-LINEAR LOSSES IN FERRITES 


I. INTRODUCTION 


The determination of the rf field strength at which nonlirear effects 


(1) 


developed the general theoretical approach and applied it to the configuration 


commence in a ferrite has been the object of several investigations. Suhl 


in which an rf driving field is applied perpendicular to a dc magnetic biasing 


(2) 


the results for a rf driving field which is parallel to and small compared to 


field. Schi@mann and his associates in various publications have presented 
the dc magnetic biasing field. This investigation considerec the special 

ease of a ferrite toroid biased to saturation by means of a pulse of current 

in a wire coincident with the axis of the toroid. Turing the time interval 

of the current puls:, the magnetization throughout the sample is circumferential 
and has the saturation magnitude, M,- Wher the current is removed, it is 
assumed that domains are formed. In eacn domain th: sigs: igation is uni- 
directional, but not necessarily circumferen’..7'. © «4S the saturation 
Magnitude. «© .@ magnetizuvion fi ~. domain to domui:, 1s not unidirectional 
however, and the average circumferential magn2tization throughout the toroid 


is reduced to the remnant magnitude, Mp: 


After the biasing cvrrent is removed, & f driving, field h, cosw, t is 
applied in the circumferential direction. The purpose of this analysis is to 
determine the minimum value of 4 which will cause instability in the toroid. 
The critical field will be a minimum in @ domain in which the dc part of the 
magnetization is in the same direction as the applied field, so that only this 
case will be considered. It is assumed that there exists at least one domain 


within the toroid for which this applies. 


(1) H.-Sulil, (Proc. DRE, 44,2270 (1956) 


(2) E. Schlomann and R. I. Joseph, Jour. Appl. Phys., 32, 1006 (1961) 
for example. 


II. DEFINITION OF TERMS AND EQUATIONS OF MOTION 


Consider a domain in which the dc part of the magnetization is in the 
same direction as the applied rf field. Assume that the domain is small 
enough to neglect the curvature of the fields in the circumferential direction 
of the toroid and let the direction of the applied field be defined as the Z 


axis of a cartsian coordinate system. 


The total magnetization within the domain of interest is assumed to be 


made up of several components defined as follows: 


is the total mangetization vector at a point 


is the dc component of the magnetization vector 


Us yxy 


is the spacial average of the ac part of the magnetization averaged 
throughout the sample 


is the iocal deviation of the magnetization vector from the spacial 
average. 


°F 


Similar definitions wili apply to the magnetic field vector terms. The 


magnetization and field vectors are then expressed in rationalized MKS units 


by 
Yi 


B 
o- x+h+m-= - (2) 
ay ~ lu ~~ Me 


(1) 


i] 
t= 

+ 
Pa=| 

+ 
te 


The magnetization is assumed to have a dc component only in the z direction. 
After the initial biasing magnetic field is removed, it is assumed that no 
applied de component of magnetic field exists. Let x and Zz be unit 
vectors in the cartsian coordinate directions and Equation (1) can be written 


in component form as 
WW = m+ im) gt (y+ om) y+ (HL +m, + fm.) 2 (3) 


The magnetic field inside the domain is made up of four terms, the 
applied rf field in the z direction, the demagnetizing field resulting from 


the average ac magnetization throughout the domain, the demangetizing 


ih 
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le gl 


NRA 
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= 
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= 
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field due to the magnetization in the toroid sample -utside the domain of 


interest, and the demagnetizing field due to local deviations of the 
magnetization. The demangetizing field due to the average ac magnetization 
in the domain can be expressed in terms of the demagnetizing factors 
NV; Bhs and N° The local deviations of the magnetization car be expressed 
in terms of a spin wave spectrum which gives rise to a magnetic field which 


must satisfy Maxwell's equations. The resultant demagnetizing field can be 


expressed ae) 
( 80-) 
hes 3 (4) 


where k is the spin wave vector. 


It will be assumed that the ecnly important local deviations of the 
magnetization are those which are driven in some way be the applied rf field 
and that these are restricted, in any particular examplc, to a narrow regicn 
of wave vector k. Therefore, Equation (4) can be replaced by the single 


term in k which dominates. That is 


k(6m-k) 
aed 


The demagneeizing ficld due to the toroid magnetization outside the 


(7 Ip 


lz-r'| 


) 
sample carn be expressed oy ‘ 


dor! eo 


~ 


hele oS) 7 


wi ww 


where ha(z) represents.the demagnetizing field at the ,oint r inside the 
sample, V represnets the gradient operator applied at th: point r and 
v' MW represents the divergence of Yfpat the point pas outside the domain of 


(3) Ronald F. Soshoo, Theory and Application of Ferrites; Prentice-Hall, Inc., 
1960, p. 232. 


(4) Ernst Schlémann; Advances in Quantum Electronics; Columbia Univ. 
Press, 1961, p. 450. 


interest. Since the magnetization is assumed to be uniform in any domain, 


Equation (6) wouid become a surface integral over all the domain boundaries. 


Since 
Paces Zs 
Vie pes) are (7) 
JES | ) [r-r'| 3 


the integrand in Equation (6) will decrease as one over the square of the 
distance from the point at which the field is being calculated to the point 
at which the source of the field is located. Thus, only the magnetization 
in domains bordering the domain in which .nstability occurs are important 

in determining h, (x). Obviously h(x) will be a function of position within 


the domain but also of the location o* the domain within the toroidal 


sample. Since the domain structure throaghout the sample cannot be known 


accurately, it will only be possible to approximate A ). The simplest 


(r 
Q ae 
approximation is to let it be zero. This is not a bad approximation for 


peints well within the domain of interest (the domain in which instability 


e 
y 


occurs) because of the one over distance squared dependence mentioned 
above. A better approximation might be to assume that the magnetization 
in neighboring domains is nearly z directed and assure a uniform value of 
by) in the z divection of the form -NM,- For most domains this would be 
a good ap>roximation since the presnece of residual magnetization implies 
an average non zero z component of magnetization throvghovt the toroidal 
sample. There will be domains, however, in which the magnetization is uct 
nearly 7 directed, and for these the approximation would be poor. In this 


analysis, it will be assumed that h,(r) is zero. 
wi 


The magnetic field in the domain can then be expressed as 


oreo UA RN a — 


K k 
H = - - = ot \ x - - ae e 
a ( Nm 2 Sm) x + ( eae = 6m-i) sy (8 ) 
k_ 
{ ’ = = -_— 3 
: sh, cosu:, t aloe fours 2 6m -k ) Z 


if, 
f 


The equation of motion for the totai magnetization vertor contains three 


terms of interest. These are: the Zeeman term which describes the coupling 


between the magnetic field and magnetization, the loss term and the exchange 
term which describes the m.tual coupling between spins. Crystalline 
anisotropy terms are umitted because we are dealing with a polycrystalline 


material. The equation of motion is given as 


a 
uae ‘ 2 es rays 2 2: ‘ 
Mi y (Ax W) i, a Hy x vat loss term (9) 
whe~: y is the gyromagnetic ratio, ue the lattice constant and Hy is the 
effective exchange field. The loss term will te of the Landau-Lifshitz 


form given by sun, 6) 


The loss term is not stated explicitly in Equation 
(9 ) because the loss parameter @ will not necessarily be the same for m 
and for 6m. Let a, represent the loss term associated with the uniform 
precession terms m and om represent the spin wAve logs term. Some care 

is necessary in the choice of sign applied to a, and @, to insure that the 


k 
homogeneous solutions (h, = 0) result in decaying waves. The following 
assumpticn is made regarding the order of magnitude of the magnitization 
M is a zero orfer quantity 
m and sy are first urder quantities 
mo bm, and ae are second order quantities 


bm. is a third.order quantity 


Since it has been assumed that the domain is saturated, the magnetization 
terms can be related by 
ye 2 


=M ( 10) 


(m+ 6m)“ + (m_ + 6m ne + (M +m + &% 
x x y y Zz 2 Z s 


Equating terms of equal order of magnitude in Equatis (10) gives tie 


approximations 
M, ~ M. 
x es 2 
Pe se 
hi 2 = 3 ~~. 
z M 
s (1) 
Sm & - ee (m 6m + m &m } 
FS M_ xX any): 


So 


ST SS eS PTR sett 


(5) Ae Gundy “Prec ERE, 4 61956) p.. 1271. 


Define w andw by 
s m 


w 8B yh 
s vor 
wo, m yM & M +m) oe yu 
Moe Yt. y(M, ye cS 
EngG: “be the 2 B70 3 “components of che fetes 


Let 7 
x 


Then from Equations (3 ) and (8 ) and the order of magnitudes assumed, 


Zand Z can be written to third order as 
or v 


. -| w,cosnyt Seu Nene Ny) ly 


Fae 
ey ma 
- E 5 |, - [wjcow,* - Wd + Wy 7 | om. 
Ke : k 
k kK ; ke k kK 
~ [Wy “2 6m +1 ¥ 2 mi. + Ty oe ]=,2, 


4 = aC Ce = N -N) 
le, osw, Wy =) m 


KK] ty 
oe = ee + [econ - we, + Wy 2 bm, 
k 
kk, Kk kk, 
+ E ‘e $m, -|Y 2 Jo, - ly Ae ] a 


(12) 


(13) 


texm in Equation (9 ). 


First Order 


Second Order (14) 


Third Order 


ES LCR RT Hu hi 


irerervtclstittstie 


First Order 


Third Order 


Second Order (15) | 
i 


Similarly the components L and a of the lcss term and E. and a of the 


exchenge term can be written 


a 
° © fe) e 
L,- ae - cht 4 M ke 
Peyisee 
1 oaks + -_— 
y™ em a, 6m, M mm 


( 16) 


(17) 


where the first term on the right is the first order term, the second tne 


seccnd order term and the third the third order term and 


2 e222 2 ,2 
: seo § 
E. a wv om, a we 
a. oe 
Ey = -@ wev bin, =a wok 6m, 


where 


wp = YH, 


(18) 


(19) 


(29) 


Tre exchange terms E. and EY are both of second order. The exchange term 


contains no first or third order terms. 


III. FIRST ORDER SOLUTION 


Collecting only the first order terms from Equations (14) through (19) 


gives 
. + = + 
(1 + ia) + iAm + ipm = iw, (cosw,t )m 
(21) 
ee a ee - 
(1 - ia )m - iAm - im = iw (cosw, t )m 
where the terms are defined as follows 
+ ‘ 
m- =m + im 
xm Y 
= 2 en i 
A= w(N, - aN. au) (22) 
B= 


z : 
2 wy (N N) 


The coupled equations of motion of Equation (21) can be decoupled in the 
homegeneous case by making the change of variables 


+ + - 


m = an - bn 
(23) 
m= - b*n + an” 
where 
cosh x = Afwo 
sinh x = B/w 
@ = cosh x/2 (oh) 


b= (1 - ia) sinh x/2 
*= (1+ ia) sinh x/2 


a. «2 2 2 
w, =A -B = wy (N, N)(N, N,) 


Neglecting all terms containing i the decoupled eruations of motion become 


n+ (4 + Aas Ge t)n* G 2B w.t)n- 
iw + @, = (4 > Acosw,t)n i a, cosw,t)n 
(25) 
vie = Os = ”s + 
n - (iw, - a A)n 2. hi ag Poona + (i >, Beosw,t)n 


a. 
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i i g 


ae 
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The hemogen2ous solution is obtained by equating the right sides of 
Equations (26) to zero ana is given by 
+ + - (iw + a Ajt 
=n e fe) fe) 
peer (26) 
- - (iw + @A)t 
e*"o fe) 


Mo =i 
3° 


The particular solution is obtained by differentiating Equations (26) 
with respect tc time, with A and ng considered time variable, and sub- 
stituting into Equations (25). The driving terms on the right sides of 
the resulting equations will cc.itain time variations of the form w, tw, 


_ 


and - 0, + w: Oniy those terms which have time variations nearly equai to 


the natural frequencies, + ws will be effective. Such terms exist only for 
ws aw, (27) 


Assuming Equation (27) and neglecting off resonant driving terms gives 


= ae Be pie act - 2p )t 
ane ~(2)(B) nv e HM = PuQ)t 
2 (28) 
ee w : een 
ns eee ‘i ei, ao 
fe) w. 2’ 0 


Differentiating Equations (28) again with respect to time and rearranging 


terms gives 


Rast i(w, - 2w )N ae @) (-8)n a) (29) 
oo a: ORO 2 wee 
Assuming a sciution of the form 
eee eee (30) 


where w is an cscillatory part and B a damping part of the time variation 
and substituting into Equation (29), gives, by equating real and imaginary 


terms, 


(31) 


The two solutions for n -. and for ee then consist of a low frequency oscil- 2 
lation with an exponentially increasing amplitude for one and an exponentially Py 
decreasing: amplitude for the other. = 
From Equation (25) it is seen that instability will occur if n= . 
grows more rapidly than the exponential decay terms in ree The threshold = 


of instability is given by 
2 2 
Bo = (@ A) (22) 


Substituting the value of gc from Equation (31) and the value of B form 


Equation (21) into Equation (32) and solving for w. gives 


Ww 
16(>>) 2 , 
ee a oie ow) aes ; (33) 
oe (ae ae) (5 es oa( 5 (Nea) 
vis ee 


U 


“ia 


The minimum instability threshold occurs when w, = ew: For this case we 


have 
3 
1) 2 
: 1 - (4)  - 0) 
w y bq 
iy 7 : | 
Om = 1 (34) 
H ; M 
erit (—)(N_ - N_) 
Wie x 
where ow, 
4H = 
Y (35) | 
is the low power resonance linewidth, 
Let N = O and recall that 
Nes Nee (36) 


w 
Then solving for (Ny - N) in terms of (4), using the definition of 0 


minimum critical field as 


10 


from Equation (24) and the assumption that w, = aw, gives the normalized 
3 
E 


5-342 (>) + 4) 


‘ ss 8 
a 6 oo 
‘erit %, w, is 37 
Gy Ue = 1 2 oo 


The critical field obtained ty this approximation is the field which 
will cause the entire domain to go unstable simultaneously. In this 
approximation, the spin system throughout the domain is in synchronism. 
For this model the critical field approaches infinity as the ratio () 


1 
approaches 0.5 and is imaginary for smaller values of (wy Ww Since finite 


ay 
critical fields have been measured below (w,,/w, ) = ea, ae model is 
obviously too crude to be useful in the range of small (w,,/,)- As the 
ratio (w/w, ) becomes larger the critical field given by Equation (37) 
grows smaller and the approximation is somewhat better. A more accurate 
approximation of the critical field for all ranges of (w,,/, ) is given by 
the second order differential equations which are taken up in the next 


section. 


“ 


a 


SS 

IV. SECUND ORDER SOLUTION = 

Retaining the definitions of Squation (22) and including first and = 
second order terms from Equations (14) through (19) gives the coupled ze 


equations of motion as 


a 


. = + =_ + — 
(1 + ig )mm + Cure ia, ) 6m— = + (iw ,cosw, tm + iAm + iBn’ 


(38) 
+ 7, a Soni Bin ae etre 
(iw cosw, t ) dm + iA + iB — bm 
where 

Bi geawh MMNe kere pia yn Phi eee aue aecee 4 
A, - dy + 28 wk ions 2 = - wy, +a wk +35 sin 8 : 
w w 3 
+ M + Dine er 2 2 
Boe (xt) = et AP ginte (39) ' 


The angle @ is the angle between the z axis of the coordinate system. and 

the direction of the spin wave vector kK. The angle » is the angle between the 
projection of k on the xy plane and the x axis. Equation (38) can be solved 
approximately by applying the technique used in solving the differertial 
equations in Part II. First the homogeneous equations for mr and om are 
solved separately. These solutions are then inserted into the particular 
equation with the constants from the homogeneous solution allowed to be time 


variable. Next all driving terms which do net have time variations nearly 


equal to the resonant frequencies from tiie homogeneous solution are 
The resultant equations can then be decoupled by a transformation 


discarded. 
similar to the one in Part II and the resulting equations solved directly. 


+ : : : A 
The homogeneous solution for m— was obtained in Part II and is given 


+ 
py Equation (26). The homogeneous equations for 6m- are given by 


(1 ¥ ia )on5 = + UA, bo + SB bm (y@) 
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Making the transformation of variables 


+ + = 
6m = ay, 6h = db, $n 
(41) 
6m = -b, 6n +a, én” 
Teese reve ara “pangs 
gives 
eS eee (42) 
6n- + (iw, + a, A.) $n- = 0 
where 
cosh x, = A /w, 
sinh x, = B /w,. 
(43) 


a, = cosh x, /2 


(1 + io) et * sinh x,/2 


low 
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Beaty Gora ae Wr ates 
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Equetion (yo) has solutions of the form 


+ + +(ib + 
6n- = én— et (li, oA )t (4h ) 
+ = 
From Equations (23) and (26) it can be seen that m andm have frequency 
components at +w and at -w . Similarly, from Equations (41) and (44), 


+ 
6m— have frequency components at tw, and at -w Assuring that these terms 


k a 
dominate, consider the frequency components of the terms in Equation (38). 
These are summarized in Table I. The equations of motion are now best 
solved in approximate form by separtely considering various cases in the 


frequency relationships. 
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TABLE I. APPROXIMATE FREQUENCIES OF TERMS IN 
SECOND ORDER MAGNETIZATION EQUATION 


A. Case 1: w, # ws dw, SF ww, fw 


From Equation (2h) it is seen that w is a function of the ferrite 
material through Wy and a function of the domain shape through the demag- 
netizing factors NT, Le and No. For a given domain then. Wo is a constant. 
The driving frequency, W,. is controlled so that the first inequality of 
this case can always be established. Spin waves in the sample can exist 
over a very broad frequency spectrum. The approxiuate equality w= 30, 
will always be satisfied for a portion of the spin wave spectrum regardless 
of the value of w) and will be the dominant term. Other spin waves with 
frequencies Ww) + Wy will not be driven by the epplied rf field and can 


be neglected. 


Considering only those terms of Equation (38) which vary at w. as 


determined from Table I gives 


(1+ ia, ) b= = + iw cosw, tém + iA bm a an (45) 
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Again making the change of variables indicated in Equation (41) and 
g 


: : ; ; ase : Be : 
assuming a solution of the form given in Equation (44) with Sn time varying, 


gives 
+ a3 s i209 i(w, - aw, )t 
6n_ & - 5 \5,| is e a k 6n 
(46) 
now 5 IRL ibe - (w - Jt 5, + 
Ww 
k 
Taking second time derivatives and combining equations gives 
2 {3 | w 2 
3 ts k S. + 
ee Rr ee ae (47) 
we aaa), 2 wy o 
Assuming a solution of the form 
én * 2 on * ot iiw + K)t (48) 


where w is the oscillatory and X the loss parts of the time variation and 


solving for w and X gives 


w= (-4,) 
Bl] w 2 w 2 (ug) 
K* = (ad Sees eae w,) 


The instability threshold for this case occurs when 


2) 


2 2 2 (50) 
KX” = a A, 

: : , Caos 
Solving Equation (50) for w. gives 
2 
2 k aX, Seal 
) 2 ( a= w,) i ey a 


el [2] 


a (51) 


Sones) 
w. crit = (yh, 
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The minimum critical field occurs when w 


AD. . 
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(ea aes 
Nhe he es Bin’ 6 
where 
2a, w 
kk 
ae 
aH = V 


is the spin wave linewidth. Using the definition of Wy 


i 


= eu, which gives 


(52) 


(53) 


from Equation (43) 


2 2 
i = 2 = 3 + 
and the assumption that w, = 2a, (a wk wy) can be eva’ uated and 
substituting into Equation (42) +o obtain the result 


GR) ee : YM, 


a 


a 


This expression is plotted in Figure l. 


( 54) 


( 55) 
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Case 2: w) - Ww = uw, 4 w 


The equations oft motion for this case are given by 


a + a, (56) 


‘+=, eee, ee eee + 
ey) 6m— + iA 6m- + 1B — bm = + iw cosw,t m- 


which correspond to a coupl.d pair of lossy harmonic cscillators on the 

left with a harmonic driving term on the right. The equations can be 
decoupled by the coordinate transformation uf Equation (41), however. one 
need rot do so to determine that She colutions of Equations (56) are 

statle ht~monic oscillators as long as the right sides of the equations have 
constant amplitude, i.e, mr constant in amplitude. The critical field for 
this case «ould then be given by Equatio: (37) which is the Field required 


te make mm grow exponentially. 


rope 
€ 


Case 3: Wi = Wo = 5 1 


For this case all of the terms from Equation (38) must be retained in 
part ard a direct solution is too comple. algabraically to be solved. A 
reaso:.able approximation to this case might be given by Equation (55) 
for low values of Wy MD and by the lower of the values from Equations (37) 


) . £ 
and (55) for large values of Wyy/ Ws - 
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V. DISCUSSION OF RESULTS AND HIGHER ORDER TERMS 


The first and second order solutions both indicated that the normalized 
critical field increases monotonically as the ratio w/o, decreases. The 
first order solution does not have a finite critical field for wy /Wy “ess 
than 0.5 and is thus inadequate for low values of Wy/, - The second order 
solution has finite values for the normalized dritical field fcr all non- 
zero values of Wy/, « For this reason, and because it considers more of the 


interacting energy terms, the second order solution is to be preferred. 


There are, however, several possibi! ‘ties for error in the second order 
solution. One of the most serious problems is the normalizing term An which 
is the spin wave linewidth. This term is related to the spin wave loss or 
relaxation, and is expected to be a function of frequency and material and 
possibly even the shape of the aample and domain rather than a constant. 
Until this term is known more explicitly, Equation (55) and Figure 1 can 


only be considered as rough approximations. 


It was assumed that the demagnetizing field due to spacial average 
magnetization in the domain couid be expressed in terms of the demagnetizirg 
factors Ny a and No These factors apply precisely only for ellipsoidal 
samples which are small compared to the wavelength of the externally applied 
fields. For other shapes constant N factors can be assumed throughout the 
bulk of the domain which will closely approximate the physical situation. 
However, near the domain boundaries and particularly near sharp edges, the 
demagnetizing field can be a rapidly varying function of position. This 
can be accounted for by allowing the demagnetizing factors to be functions 
of position near the domain boundaries and constant throughout most of the 
qomain. Unfortunatel:. computation of the functional dependence of the 
demagnetizing field (or factors) is a prohibitively difficult task for 


almost all shapes of interest. 


The effect of th2 v riation in N, can be seen from the spin wave dis- 
persion diagram shown in Figure 2. The abscissa scale in Figure 2 is 


{ > 
normalized by the factor if w,4/ 8 w, which depends on the material constants 


E 
Vos Me: a, and He and which makes Figure 2 a general curve applicable to 


all ferrites. 
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For a fixed value of w/w, (recalling that the applied Ww) = au) it is 
seen that for various values of N, spin waves of various wave numbers or wave 
lengths can be generated in di*ferent regions of the sample. If N, varies 
rapidly enough to give interaction between the local demagnetizing fields, 
6h, Equation (5) may no longer be valid. 1n this case Equation (4) would 
not be valid either since it represents the case in which a spectrum of 
spin waves is generated in a given local region. For the case in which 
spin waves of differcit wave numbers k are generated in neighboring regions 
within the domain Equation (4) should be replaced by a weighted sum in which 
the weighting terme would be at least a function of the distance between 
regions. It can also be noted from Figure 2 that for small values of uy, /udy 
and for some angles g two spin waves with the same frequency can he 
generated with different values of wave number k. This is possible for 
values of g which satisfy 


sin’p <N, (57) 


w 
and values of ze which satisfy 


Eo (58 ) 


This situation again makes Equation (5) invaliu unless a second term is added 
to account for the degeneracy. This situation is not too serious however, 


since it is expected ty /oy > No for most practical cases of interest here. 


So far only the first and second order solutions to the instability 
threshold problem have been considered. Some qualitative statements re- 
garding the third order solution can be made. The exchange terms contain 


no third order quantities and so remain unchanged. The third order less 


eb 


terms are off frequency and can be neglected. The third order Zeeman terms 
allow direct coupling from the average magnetization terms to the local 
deviation terms. Neglecting off frequency terms, the third order equation 


for the magnetization can be written as 


se . k 
(1 - ia, ) sm = 1A, gut 1BY gm =i fescose,t Bocatata")| én* 
k 


- iy 4S om (59) 


If the driving term containing W, is much larger than the remaining driving 
terms, Equation (59) reduces to the second order solution. If the term 

w. nesligibly small compared to the other driving terms, Equscion (59) will 
reduce to Equation (21) of Suhl's analysis. In general, the local spin 
system can gain energy from the driving field through the W term or gain 
energy from or lose energy to the spacial average spin system through the 
cross product terms. A complete solution should indicate that each driving 
term will result in one exponentially increasing local spin wave and one 
exponentially decreasing one above some threshold value. If the wave which 
grows from the we pumping term corresponds to the decaying wave from the 
cross product terms, the threshold would be e.-pected to be higher than 

that computed in the second order solution resulting from the transfer of 
pump energy from the local spin deviations to the spacial average. This 
would tend to distribute th. total spin energy in the sample more uniformly 
and thus seem to be a correct evaluation from a thermodynamic approach. 
Except in the case where w) + eu, as well as Oe eu, an ancrease in the 
spacial average magnetization due to direct pumping is not expected. In 
this case a noticatle increase in the average magnetizatior due to indirect 
coupling from locally deviating spin system which is being pumped would not 
be expected because of the large difference in the numbec of spins involved 
in the local variation as compared to the domain as a whole. Therefore, 
except in the special case of pumping et twice the bulk i1esonance frequer cy 
of the domain Ww,» the normalized critical field can be expected to be scne- 
what larger than that computed in the second order solution due to the 
third order Zeeman terms and the instability threshold will be caused by the 
direct coupling from the applied rf magnetic field to the local spin system. 
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PART III. FROPAGATICN OF TE MODES IN DIELECTRICALLY LOADED WAVEGUIDES 
ABSTRACT 


The propagation of TE no modes in rectangular waveguides which contaii 
two dielectric slabs parallel to the narrow wall and extending over the full 
height of the guide is investigated. Waveguide and dielectric are assumed 
to be lossless and infinitely long. Apart from these restrictions the 
dielectric slabs may have arbitrary thickness, position and dielectric 
constant. The analysis is restricted to TE no modes with the E-field parallel 
to the narrow guidewall. The guide containing only one dielectric slab is 
covered by this analysis. The even modes n = 2,4,6, of the guide with two 
slabs correspond to the odd modes n' = n/2 = 1,2,3,... of the guide with one 
slab and half the width of the guide with two slabs. 

For the relative dielectric constants € = 2.25,4,9, 12.25, 16425 the 
normalized chtoff frequencies for TE 10, 20, 30, 40, 60 modes have been 
computed for fifteen slab thicknesses (including O and 100% filling factor) 
and a maximum of eleven slab positions, normalized propagation constants have 
been computed for five slab thicknesses and a maximum of eight slab positions 
for TE 10 and TE 20 modes between their respective cutoff frequencies and 
a frequency slightly above the second and fourth order mode cutoff fre- 
quency of the Empty guide, respectively. The results are presented graphically 
and numerically. 

These results &re discussed. The parametric dependence of field dis- 
tributions, the ratio of magnetic field components (ellipticity), the 
ratio of cutoff frequencies (fractional bandwidth), and normalized wave 


impedances are illustrated. 
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PROPAGATION OF TE MODES IN DIELECTRICALLY LOADED WAVEGUIDES : 
INTRODUCTION { 

Previous analyses of propagation of TE no modes in rectangular waveguides 
which contain dielectric slabs have in general been restricted to two cases, pe 
where the dielectric slab is placed a) against a waveguide wall (1, 2, 3, 4), 4 
or b) in the center plane of the waveguide (2,3,4,5,6). A few more special 4 
cases have been considered by workers dealing with ferrite applications in . 
the microwave region (4, 8). For a rather general position of the dielectric . 
slab certain phase shift characteristics of the loaded guide have been : 
calculated (7). ; 
The present analysis deals with a rectangular waveguide which contains : 

F 


two dielectric slabs parallel to the narrow walls and extending over the full 
height of the guide. The slabs are placed symmetrically with respect to the 
center E-plane of the guide. Apart from this restriction the slabs have 
arbitrary position, thickness and dielectric constant. Only TE -modes with 
E-fields parallel to the narrow guide wall are considered. The guide and the 
diele-’ -ic are assumed to be lossless and infinitely long. Waveguides 


containing only one dielectric slab are covered by this analysis. 


QUANTITIES AND SYMBOLS 
Figure 1 shows the crods-section of a rectangular, dielectrically loaded 
The broad dimension of the guide 


waveguide in a rectangular coordinate system. 
extends along the x-axis, y is the direction of propagation of fields in the 
guide. a, c, d, w are waveguide dimensions along the x-axis. B = w Vt5%o 
is the free-space propagation constant, k is the propagation constans in the 


ie 


wilt 


guide in the direction of the guide, p in the empty guide region and q in the 
dielectric are propagation constants transverse to the direction of the guide 
Instead of these symbols dimensionless quantities 
These are obtained by either multiplying 


and tne electric field. 


will ve used throughout the analysis. 


or dividing the above given quantities by w. Quantities, their dimensions 


2 

z 
i 
| 
5 
3 


(MKSA-sys and symbols used here are: 
iE 
: 
8 
h 


Symbol Unit 
a = a/w S 
6 = d/w = 
Y = c/w = 
pea a m 
p= x/w : 
Loiyet, =; 
t Ss 
-1 
Ww Ss 
B = Bw - 
B = 6 w - 
c c 
K = kw - 
Ps pw = 
Q = qw = 
Py Py radian 
6 radian 
C,D - 
\ vsA~4m7+ 
. Seek 
€ AsV “nm 
fe) 
€ or DK - 
B=i.E vm> 
Eo Vas 
Hei H+i # ne 
XOX oF 
He AM 
Y avi 
u -1 
XY AV 
re) 
* - 
ELL - 


(y) Raw d = ce Cee d—wea-= = (x) 


Quantity 


waveguide dimensions as shown in Figure 1 


a+6+ye=1 
right hand coordinate system as shown in Pig. 1 
normalized x-coordinate 
unit vectors 
time 
angular frequency 
free-space propagation number (frequency parameter) 
normalized cutoff frequency 
longitudinal propagation number in the guide 


transverse propagation number in empty part of 
guide 


transverse propagation number in the dielectric 
electrical widths of waveguide sections 

phase angle 

relative amplitudes 

free-space permeability 

free-space dielectric constant 


relative dielectric constant of dielectric 
(DK appears in the computer printed tables) 


electric Pied 

normalizing electric field 
magnetic field 

normalizing magnetic field 
wave admittance of loaded guide 
wave admittance of free space 


fractional bandwidth 


ellipticity 
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Figure 1. Waveguide with Dielectric Slabs 
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THEORY 2 

It completely suffices, because of the symmetry of the loaded guide, to = 
consider the regions: 1, 12, Lift of Figure 1, so that 0 = x 5 woor O°s' 9 Sil. a 
Atl fields vary as exp fwt - ky), so that 0/dt = jw and 3/dy = -jk. This +t 


t a ty dependence is omitted in ail equations. The relative permeability 
of the dielectric is assumed to be 1. Maxwell's equations for the problems 


considered here reduce to 


Be Be esi a nas ; : : 
(1) oe Spirty ag odie, (iH + iH) 
oH oh, 
(2) D7 re ye 
(3) VE = 0 
(4) VH = 0 
and from equation (1, 2, 3) 
; Ae 
(5) VvxVx B= - (25+ 25) b= wees 
ox oy 

The boundary conditions are 
(6) E = E is continuous, H = H_ is continuous 

tang tang y 


The fields in the various regions of the guide can be described as in Table 


1 with dimensionless quantities: 


» E/E, JH ek = Ha /KE,, 
Posk/pst | K/B-1 | 1<K/as/@" ae £ 1<K/B:/e 


Cch |P| (1-p) re 0 -C Fn [Pl (1-9) 


1-Y<)sl1_ | even  coinP(a-9) C(1~p) | Csh iP| (1-9) -CEcosP(1-9) CR oF en ipl (1-9) 
Table 1. Normalized Field Distribution in Waveguide which Contains Dielectric Slabs 4 
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2K/3 = 3.33 


2K/3 = 3.14 


Normalized Field Distribution in Waveguide with Dielectric Slabs 


Figure 2. 
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These field distributions are illustrated by Figure 2 for the first odd 
(n = 1) and first even (n = 2) mode. The even modes (n = 2,4,6 ....) of 
the waveguide under consideration with two symmetrically placed dielectric 
slabs, where 2w is the waveguide width, correspond to the odd and even modes 
(nt = n/2 = 1, 2, 3, ....) of a waveguide with only one dielectric slab, 
where w is the waveguide width; that is the left half of Figure 1 only. 

The fields of Table 1, together with equation (5) give for the propa- 


gation constants the relations 


5) a eae ere 
oe 


(7) Pe us Ke 
(8) Q sen oe 
(9) Poce co mea(er a ens 


For any Qe > O there exists a pe 2 O, the sign depending on (e - i ae 
Qe > O describes a sinusoidal field distribution in the dielectric slab, 
pe > 0 describes a sinusoidal field distribution in the empty waveguide 
region. For ce < 0 the field distributions in the empty waveguide regions are 
hyperbolic functions which describe a quasi-exponential decay. p = 0 gives 
the intermediate function where E and HY have a constant slope in the empty 
part cf the guide ard where we is constant. In this particular case the 
fielus in the empty region are described by functions as » exp(Rp) = 9, 
since P= 0. When QF <0, also P <6, Then tne boundary conditions (6) 
are violated. Hee is no longer continuous. Therefore, Qe 2 0 for all TED 


modes. oe = 0 is reached at infinitely high frequencies. The assumptions of 


Tacle 1 cover cnly and all allowed field distributions for TE 6 modes and 
are thus justified. While the frequency increases from cutoff to infinity, 


K/B increases from 0 to/€. For K = B the field distribution between the 


i 


two dielectric slabs represents a pure TM field for all odd TE Ao modes in 


the waveguide considered. 


The determinental equations for B{K) - e.g. cutoff frequencies for 
k = O - or K(B) - propagation constants for given frequencies - are found by 
expressing the widths of the empty sections of the guide in equivalent widths 
of a guide completely filled with the dielectric. The field distribution in 


the slabs is not changed by this replacement. The total electrical wilth 


TRAMP. amanactor Hannan neh 2LNt HanHT HRCI aig 


(2w,) of this equivalent guide is nn for a Teo mode, 


« 


:3 
r= 
= 
= 


Av the boundary between regions I and II of Figure 1, Table 1 and 


equation (6) yield for frequencies © s K/R< 1 and with p, = Fy, P= ga + 8 


lu) E si % sin 
( ! ae in P, E. Dr sina 
C8) jELP ces P/M. = JE. DQ cos p./We 
i \ far ie) Co 
c ) od H ih \ = jee oe +f =- j —— neyse 
(12) 2, (a B(o)/JK, (a) = 54> te ey = J am te 5 


a /P = a5 and Oh /s Zo are the transverse wave-impedances of the empt.s 
and loaded waveguide regions, respectively, for waves travellinyg in + x 
directicn. 2, (a) is the impedance experienced by a wave travelling into 
the shorted waveguide of impedance Zp and length p, or impedance an nd 


wength p,. From equation (12) follows the equivalent length 


\ mone peat ae’ 
(13) p, = arcte(S te p,) 
Similar considerations for higher frequencies, odd and even modes, and hot 
slab boundaries lead for the various waveguide regions to the actual and 


equivalent electrical widths of Table 2. 
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Widgy Be Queney | OSK/B<1 K/B= 1 1<K/Bs/€_ a 
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Equivalent I p. | arcte(# tg p,) : arctg(Qp, ) arcte(, tho, ) 
II P| 26 es ee: 

; III PY ride 7 : PLY. 

eee ee eee ee : odd modes even ee aa ame fey 

mene | ers a) ii 

jEquivalent III p. arecte (B) tap, |: acrte(Qp, ) arctg(-1) 7 tho,, | 
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Table 2, wlectrical Widths in Waveguides, which Contain Dielectric Slabs 
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4 
Pe is defined for 0 s K/B < 1 for odd modes by 1/Ptg P= 1/Qte P. and for even : 
modes by tg p/P = tg a./Q and for other frequencies accordingly. Table 2 is s 
the skeleton for the computer program used to determine the cutoff frequencies 2 
B(K = 0) and propagation constants K(B). The determinental equaiton is in 4 


both cases 


TT ABIMI ta 


y = % 
(14) Pe = nn/2 for oe -modes 


Applying equation (6) to the boundaries of regions I/II and II/III of 


Figure 1 yields the phase angle 8 and the relative amplitudes of Table 1. 


They are given in Table 3 in terms of the electrical widths defined in Table 2. x 

= 

ee Sse are se a 

sin(p, + P,) sin(p, + p)/ch/o), | 2 

| eS 

sin(p, + p,)/sinp,i sin(p, + 3)/P | sin{p, + p)/shp,, | . 

Table 3. Phase Angle and Relative Amplitude of Fields in Waveguides, E 

Which Contain Dielectric Slabs M] 

COMPUTER RESULTS E 

y 

Normalized cutoff frequencies have been calculated for TEN modes of : 

the guide of Figure i with n= 1, 2, 3, 4, 6 for six relative dielectric E 
fonstants \e.eo, 43 9, .25, 16, 25), fifteen slab thicknesses (including 

0% and 100% filling factor) and a maximum of eleven positions of the slab 

in the guide. Normalized propagation constants have been calculated for TE) 6 : 

and TES, modes between their respective cutoff frequency and a frequency i 

somewhat above the cecond and fourth order mode cutoff frequency of the enpty = 

guide, respectively, again for the relative dielectric constants given E 


above,five slab thicknesses (5, 10, 15, 25, 40% filling factor) and a 


= 
= 
= 
= 


maximum of eight slab positions. The position pevameter a + 6/2 gives tne 
distance between the left guidewa!l and the center plane of the left slab 
as a fraction of half the guide width. For the odd modes the position is 
varied between the slabs touching the guide walls (a + 5/2 = 6/2) and the 
slabs touching each other in the center of the guide. (@ + 6/2 =1- 6/2). 
For even modes it suffices, because of the symmetry of the field distribution, 
to vary the position between the slabs touching the wall and woving them haif 
way towards each other (a + 6/2 = 0.5). For even modes the cutoff-frequencies 
and propagation constants are the same fora@ + 6/2 = 7 anda+ 6/2 = 1 -7. 
The normalized cutoff-frequencies are given graphically in Figures & ... 13 
and numerically in Tables 4 ... 9. In these tables the slab thickness (4) 
is the horizontal. parameter, the rositicn of the slab (w + 6/2) and the 
order of modes (N) are the vertical parameters. 

The normalized propagation ccnstants K are given graphically in Figures 
14 ... 43 with one set of curves for every dielectric constant and slab 
thickness combination ard the position as parameter for every curve. They 
are given numerically in Tables 10 ... 29. The horizontal parameters are 
position (@ + 6/2) and order of mode \TELG) Tne vertical parameter is the 
normalized frequency (B). The normalized cutoff frequency (2) for each 
set of parameters is given cn top of each column. Propagation constants 
-0.0 indicate that the frequency for this K is below cutoff. 
Examples: 
a) Guide WR 137, width 1.372 inches, two slabs each 0.C69 inches thick, 
0.206 inches between left wall and center of left slab, relative dielectric 
constant € = 9, w = width/2 = 0.686 inches, 6 = 0.1, a+ 6/2 = 0.3. Wanted: 
TE) 6 and TES5 cutoff-frequencies and guide wavelengths at 5.46 GHz. One finds: 
B (TE) ,) 1335 alae 2.05. With A, = arw/E one gets d ATE 
KES) 5.35 em, f (TE 
one finds B = 2 and K(TE, 


10! SME oh an 
2¢ 3.63 GHz f (rE,,) 56 Gi, At 5c.5 GHz 
) ; 


- } 3 corresponding to A_ = omw/K 5.74 em, 
& 


no propagaticn for the TE, mode. 


b) Guide WR 90, width 0.9 inches, one slab with € = 12.25, 0.135 inches 
thick, 0.18" between wall and center of slab. Wanted: cutoff frequencies of 
the two lowest order modes, guide wavelengtr at 10 GHz for the lowest order 


a (ee ees ee . ss an a Ze CUE On Loss 
mode .w' = width = 0.9 inch 3, 6 = 0.15, a + 8/2 = G.2. B (TE, ,,) a B (TE, .) 


of guide with two slabs and with 1.8 inch>s and w = width/2. Cue finds 

B (TES) = 1.83, B A (TE54) = = (TE, ,) = 4.25) witn do = 2mw'/B one gets 

ho (TE, i) = 7.85 om, A (TES im 3.38 om and (TE, . ) a 3.82 GHe, f o (TS, oo 
7.83 oh, At 10 GHz one Fait B= 2mw'! |r, = 4.8 ena Ke d2ena5 ies 
A =: 1.12 om. 
8 


ce) The TE, solutions with a + 6/2 = 0.5 are equivalent to the TE 6 


solutions with the two slabs touching each other, where a + 6/2 = 1 - 6/2, 
me. 2h (TE, |: B,@+ 6/2 =1 - 6/2) = K (TE Eno) 2B, & + &/2 = 0.5) 
2B, (TE, a+ 6/2=1 - 5/2) = B (TE, a4 6/2°5 0.5) 


: 
i 
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LISCUSSION OF RESULTS 


The influence of the dielectric slab on the field distribution in 


TATA 


the waveguide is te concentrate within the slab with increasing frequency 
an increasing fraction of the total energy flowing through the guide. The 
phase velocity approaches asymtotically "/@ times the velocity of light in 
free space". The paranetric dependence of the E field distribution is 
illustrated by Figure 3. 


The influence of a thin slab on t¢ie cuteff frequencies is the stronger 
tne closer the slab position is to the lines of maximum electric field 
Strength jin the empty guide. With increasing slab thickness the influence 
is less and less related to the empcy guide field distribution. For thick 


Slabs it is weakest when the slabs touch the guide walls; for odd modes 


“vn HR 


it is strongest when the slabs are in the center of the guide and touch each 


other, for even modes it is strongest when each slab is in center of half 
& guide width. 


At a certain thickness the cutoff frequency is nearly independent of 
the slab position (see TE 36, ho, 60?" For thinner slabs a similar behavior 
can be found in the propagation characteristics. Figure 4 shows K versus 
B curves for the TE, 6 mode with € = 9 and &= 0.2. At B. 5.42 the prop- 


agation constant is independent of the slab pesition. Above this crossover 
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frequency the same order of curves is obtained as for eutoff frequencies 


PARAMETER: a +8/2 = 
ot | ttt yt 


Figure 3. Propagation Characteristics in Waveguide with Dielectric Slabs 


« = 9.0, a +8/2= 0.5, PARAMETER : 8 


8=0.2,a +8/2=0.5, PARAMETER: € 


€ =12.25,8=0.1, PARAMETER: a+8/2 


E-Field Distrubition in a Waveguide with Dielectric Slabs 
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with thick slabs. The propagation characteristics for TE, modes for 

e.g. © = 12.25 and & = 0.1 show, however, that this order, too, may change 

at even higher frequencies. If the empty regions in half the waveguide 

(I and III, w and y in Figure 1) are not equal, a wave may treat tiem at 

high frequencies as essentially equal with the value of the narrower region. 
The waveguide shows then a larger, propagation constant and field concen- 
tration at a given frequency than a guide with @ = y does. An E-field distri- 


bution for such a case is shown in the lowest part of Figure 3. 


The ratio of the magnetic field components, the ellipticity, 
(15) ELL(p) = H, (@)/5u,(®) 


at the slab boundaries I/II and II/III approaches “one" asymtotically, as 


illustrated in Figure 5 for TE_. and TE 0 modes for various parameters. 


10 2 
How the ratios of cutoff frequencies, the fractional bandwidths 


(16) = Ealan0. and = BCPFo) 
\ Tea eG Sale eG) 


vary with different parameters is shown in Figure 6. Finally, tne parametric 


dependance of the normalized wave admittance 
= = = 1 = ‘ = aia 5 -1 
(17) Man (H,/E,)/¥, = K/B Ah, with Y= /€7k5 = (377 ohms) 
is shown in Figure 7. 
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